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We present the next-to-leading order calculation of two four-jet event shape variables, the D 
parameter and acoplanarity differential distributions. We find large, more than 100 % radiative 
corrections. The theoretical prediction for the D parameter is compared to L3 data obtained at the 

peak and corrected to hadron level. 
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In the second phase of the Large Electron Positron 
Collider it is an important question how well the char- 
acteristics of QCD four-jet events, i.e. events in which 
an s-channel Z*^ or 7* decays into four quark and gluon 
jets, are understood at large energies. This question is 
of interest because W~^W~ events lead to four-jet final 
states for which the main backgrounds are QCD events 
and because QCD four-jet events are also the principal 
source of background for Higgs and other new particle 
searches. The perturbative description of the QCD four- 
jet events is also interesting in its own right as a tool 
for testing perturbation theory in a regime with small 
hadronization uncertainty and for measuring the QCD 
color charges j^, or as a means of testing whether ex- 
perimental data favor or exclude the existence of light 
gluinos d). 

Recent theoretical developments make possible the 
next-to-leading order calculation of four-jet quantities. 
There are now several general methods available for the 
cancellation of infrared divergences that can be used for 
setting up a Monte Carlo evaluation of next-to- leading or- 
der partonic cross sections ■ The main ingredients of 
the calculation are the four-parton next-to-leading order 
and fi-ve-parton Born level squared matrix elements. The 
tree level amplitudes for the processes e+e~ —> qqggg and 
e"'"e~ —> qqQQg from which the latter can be constructed 
have been known for a long time Recently Camp- 
bell, Glover and Miller calculated the other vital piece 
of information, the virtual corrections for the processes 
e+e^ ^ 7* — > qqQQ and qqgg [Q. Also, the new tech- 
niques developed by Bern, Dixon and Kosower in the cal- 
culation of one-loop multiparton helicity amplitudes |^] 
made possible the derivation of analytic expressions for 
the helicity amplitudes of the e"*"e~ — > , 7* — > qqQQ 
process 1^ and results for the other subprocess are ex- 
pected to appear soon. Using these results Dixon and 
Signer calculated the next-to-leading order corrections 
for four-jet fractions with various clustering algorithms 

In this Letter we enlarge the list of four-jet observables 
that are calculated to next-to- leading order accuracy. We 
present results of the calculation of QCD radiative cor- 
rections to two four-jet shape variable differential distri- 
butions — the D parameter and acoplanarity. 

We use the matrix elements of Ref. j?) for the loop cor- 
rections. In the calculation of these matrix elements all 
quark and lepton massess are set to zero, thus our results 
are valid in the massless limit. We note that the results 
in Ref. do not include the "light-by-glue" virtual con- 
tributions which were shown to be negligible [ pT| . 

The higher order correction to the leading order par- 
tonic cross section cr'"'^ is a sum of two integrals, one of 
the real correction da^ that is an exclusive cross section 
of five partons in the final state and the other of the vir- 
tual correction da^ that is the one-loop correction to the 
process with four partons in the final state: 
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The two integrals on the right-hand side of Eq. (|l|) are 
separately divergent in d = 4 dimensions, but their sum is 
finite provided the jet function J„ defines an infrared safe 
quantity. Therefore, the separate pieces have to be regu- 
larized. We use dimensional regularization in d = 4 — 2e 
dimensions, in which case the divergences are replaced by 
double and single poles in e. We assume that ultraviolet 
renormalization of all Green functions to one-loop order 
has been carried out, so the poles are of infrared origin. 

There are several ways of exposing the cancellation of 
infrared singularities directly at the integrand level [||-D . 
The method used in the calculation presented in this Let- 
ter is a slightly modified version of the dipole formalism 
of Catani and Seymour that is based on the subtrac- 
tion method. The general idea of the subtraction method 
for writing a general-purpose Monte Carlo program is to 
use the identity 
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where dcr^ in the dipole formalism is a proper approxi- 
mation of dcr^ in the kinematically degenerate (soft and 
coUinear) region so that it has the same pointwise sin- 
gular behaviour (in d dimensions) as dcr^ itself. As a 
result, d(T^ acts as a local counterterm for dci^, that is, 
[Aa^ — Aa^] is integrable in four dimensions by defini- 
tion. The approximate cross section is constructed in 
such a way that it can be integrated analytically over 
the exactly factorized one-parton subspace leading to e 
poles, that can be combined with those in Acr^ . The e 
poles are guaranteed to cancel for infrared safe observ- 
ables (Kinoshita-Lee-Nauenberg theorem). These quan- 
tities have to be experimentally (theoretically) defined in 
such a way that their actual value is independent of the 
number of soft and coUinear hadrons (partons) produced 
in the final state. In particular, this value has to be the 
same in a given four-parton configuration and in all five- 
parton configurations that arc kinematically degenerate 
with it (i.e. that are obtained from the four-parton con- 
figuration by adding a soft parton or replacing a parton 
with a pair of coUinear partons carrying the same total 
momentum). This property can be simply restated in 
a formal way. If the function J„ gives the value of a 
certain jet observable in terms of the momenta of the n 
final-state partons, we should have 
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in any case where the five-parton and the four-parton 
configurations are kinematically degenerate. It is easy to 
prove that the observables considered in this Letter fulfill 
this property. When the requirement of infrared safety, 
relation (||) is fulfilled the second integral in Eq. (^) is also 
finite in d = 4 dimensions and a^^'^ can be easily imple- 
mented in a 'partonic Monte Carlo' program that gen- 
erates appropriately weighted partonic events with five 
final-state partons and events with four partons. 

For the precise definition of the approximate cross sec- 
tion in the dipole formalism, we refer to the original work 
of Catani and Seymour ||]. The distinct feature of this 
formalism as compared to other subtraction methods ^] 
is the exact factorization of the five-particle phase space 
into a four-particle and a one-particle phase space, and 
that the approximate cross sections provides a single and 
smooth approximation of the real cross section in all of its 
singular limits. These features lead to a well-converging 
partonic Monte Carlo program. 

In this Letter we consider two classic four-jet event 
shape variables. The D parameter fl^ is derived from 
the eigenvalues of the infrared safe momentum tensor 
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where the sum on a runs over all final state hadrons and 
is the ith component of the three-momentum pa of 
hadron a in the cm. system. The tensor 9 is normalized 
to have unit trace. In terms of the eigenvalues Xi of the 
3x3 matrix 6, the global shape parameter D is defined 
as 

D = 27AiA2A3. (7) 
The second observable is acoplanarity defined as 



A = 4] 
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where the sum runs over all particles in an event, and p^"* 
is measured perpendicular to a plane chosen to minimize 
A. 

In the case of the shape variable differential distribu- 
tions for observable O the jet function J„ is actually a 
functional, 
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where D^") is given by Eq. (0) and A^") is given by 
Eq. (I). 

Once the integrations in Eq. (Q) are carried out, the 
next-to-leading order differential cross section for the 
four-jet observable O4 takes the general form 
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In this equation (Tq denotes the Born cross section for 
the process e+e^ — > qq, /3o — ("t^^ ~ ^'^R^s) '^^^^ 



the normahzation = i in Tr(r°Tt^) = TrS"'', s is 
the total cm. energy squared, /x is the renormalization 
scale, while -B04 and C04 are scale independent functions, 
is the Born approximation and C04 is the radiative 
correction. 

The first complete results obtained for four-jet ob- 
servables at next-to-leading order accuracy are 
four-jet rates for three clustering algorithms: the 
Durham iQ, the Geneva and the EO |16| schemes 
calculated for three colors, five massless flavors and 
with as{Mz) = 0.118. We also calculated these ob- 
servables and compared the results of the two calcu- 
lations in Table I. There is a very good agreement 
between the two calculations apart from the 3 % dis- 
crepancy in the Geneva scheme result at j/cut ~ 0.05. 

TABLE I. Comparison of the four-jet fractions calculated 
by the two partonic Monte Carlo programs MENLO PARC 
and DEBRECEN (this work). 



Algorithm 


ycut 


MENLO PARC 


DEBRECEN 






0.005 


(1.04 ±0.02) 


10-^ 


(1.05 ±0.004) 


■ 10 


-1 


Durham 


0.01 


(4.70 ± 0.06) 


10-2 


(4.66 ± 0.02) 


10" 


2 




0.03 


(6.82 ±0.08) 


10-3 


(6.87 ±0.04) 


10" 
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0.02 


(2.56 ±0.06) 


10-^ 


(2.63 ±0.06) 


10" 


1 


Geneva 


0.03 


(1.71 ±0.03) 


10-1 


(1.75 ±0.03) 


10" 


1 




0.05 


(8.58 ±0.15) 


10-2 


(8.27 ±0.08) 


10" 


2 




0.005 


(3.79 ±0.08) 


10-1 


(3.88 ± 0.07) 


10" 


1 


EO 


0.01 


(1.88 ±0.03) 


10-1 


(1.92 ±0.01) 


10" 


1 




0.03 


(3.46 ±0.05) 


10^2 


(3.37 ±0.01) 


10" 


2 



We list the numerical values for Bjy , Cd in Table II and 
those for Ba and Ca in Table III. Our program generates 
four and five parton events with an appropiate weight. In 
order to obtain the Bq^ and C04 functions we calculated 
the O4 observable of each event, multiplied each weight 
by O4 and added to the appropiate bin O4. 

We define the average value of these shape variables as 
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We studied the dependence of the average value of the 
D parameter on the renormalization scale in Fig. 1. The 
strong dependence found at leading order is decreased at 
next-to-leading order. However, there still remains sub- 
stantial scale dependence showing that the uncalculated 
higher order corrections are presumably large. The fea- 
ture is similar in the case of acoplanarity, but the residual 
scale dependence is even larger. 
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TABLE II. The Born level and next-to- leading order scale 
independent functions Bo and Co- 



D parameter 



D 


Bo 




Co 




0.00 


(6.60 ± 0.02) 


• 10 


(1.08 ± 0.06) 


• 10 


0.04 


(2.32 ± 0.01) 


• 10^ 


(1.24 ± 0.02) 


• 10* 


0.08 


(1.45 ± 0.01) 


• 10 


(8.59 ± 0.12) 


■ 10^ 


0.12 


(1.03 ± 0.01) 


• 10^ 


(6.24 ± 0.12) 


■ 10^ 


0.16 


(7.74 ± 0.05) 


• 10 


(4.99 ± 0.11) 


■ 10^ 


0.20 


(5.97 ± 0.04) 


• 10 


(3.85 ± 0.06) 


■ 10^ 


0.24 


(4.69 ± 0.03) 


• 10 


(2.98 ± 0.05) 


■ 10^ 


0.28 


(3.77 ± 0.03) 


• 10 


(2.52 ± 0.05) 


• 10^ 


0.32 


(3.01 ± 0.02) 


1 nl 
• 10 


(1.94 ± 0.05) 


■ 10^ 


0.36 


(2.41 ± 0.02) 


1 nl 
• 10 


(1.59 ± 0.04) 


■ 10^ 


0.40 


(1.98 ± 0.02) 


1 nl 

• 10 


(1.37 ± 0.03) 


• 10^ 


0.44 


(1.61 ± 0.02) 


1 nl 
• 10 


(1.06 ± 0.03) 


■ 10^ 


0.48 


(1.30 ± 0.01) 


1 nl 
• 10 


(8.72 ± 0.19) 


• 10^ 


0.52 


(1.07 ± 0.01) 


1 nl 

• 10 


(7.11 ± 0.16) 


■ 10^ 


0.56 


(8.48 ± 0.10) 


1 nO 
• 10 


(5.68 ± 0.14) 


■ 10^ 


0.60 


(6.70 ± 0.09) 


1 nO 

• 10 


(4.46 ± 0.21) 


■ 10^ 


0.64 


(5.33 ± 0.08) 


1 nO 

• 10 


(3.52 ± 0.11) 


■ 10^ 


0.68 


(4.10 ± 0.07) 


1 nO 
• 10 


(2.74 ± 0.09) 


• 10^ 


0.72 


(3.11 ± 0.06) 


■ 10° 


(2.08 ± 0.08) 


• 10^ 


0.76 


(2.24 ±0.05) 


•10° 


(1.54 ± 0.06) 


• 10^ 


0.80 


(1.52 ±0.04) 


• 10° 


(1.03 ±0.04) 


■ 10^ 


0.84 


(9.95 ± 0.30) 


•10-1 


(6.66 ±0.31) 


■IQi 


0.88 


(5.74 ±0.22) 


■ 10-1 


(3.89 ± 0.20) 


■IQi 


0.92 


(2.68 ±0.15) 


•10-1 


(1.71 ±0.19) 


•lOi 


0.96 


(5.16 ±0.61) 


• 10"^ 


(2.60 ± 1.30) 


•10° 


TABLE III. 


The Born level and next-to-leading order scale 


independent functions Ba and Ca- 






A 


Ba 




Ca 




0.00 


(3.34 ±0.01) 


10^ 


(1.56 ±0.01) • 


lO'' 


0.04 


(7.39 ± 0.03) 


10^ 


(5.17 ±0.08) ■ 


10^ 


0.08 


(3.63 ±0.02) 


10^ 


(2.69 ± 0.06) ■ 


10^ 


0.12 


(2.05 ±0.01) 


10^ 


(1.56 ±0.03) ■ 


10^ 


0.16 


(1.23 ±0.01) 


10^ 


(9.59 ±0.22) ■ 


10^ 


0.20 


(7.63 ± 0.07) 


10° 


(6.12 ±0.15) ■ 


10^ 


0.24 


(4.81 ±0.05) 


10° 


(3.97 ±0.12) • 


10^ 


0.28 


(3.02 ± 0.04) 


10° 


(2.57 ±0.09) • 


10^ 


0.32 


(1.78 ±0.03) 


10° 


(1.59 ±0.08) ■ 


10^ 


0.36 


(1.08 ±0.02) 


10° 


(1.10 ±0.06) ■ 


10^ 


0.40 


(5.99 ±0.17) 


10-1 


(5.99 ± 0.33) • 


IQi 


0.44 


(3.19 ±0.12) 


10-1 


(3.67 ±0.25) ■ 


IQi 


0.48 


(1.51 ±0.09) 


10-1 


(1.90 ±0.99) • 


IQi 


0.52 


(5.91 ±0.52) 


10-^ 


(8.45 ± 0.99) • 


10° 


0.56 


(1.55 ±0.26) 


10-^ 


(2.84 ±0.42) • 


10° 


0.60 


(1.33 ±0.84) 


10-^ 


(7.64 ± 1.77) ■ 


10-1 


0.64 


(0.00 ±0.00) 


10° 


(5.55 ±2.66) ■ 


10-^ 
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FIG. 1. Renormalization scale dependence of the average 
value of the D parameter, = [ij \fa 

The same conclusion is drawn if we look at the depen- 
dence of the K factors on the observables as depicted in 
Fig. 2. In case of the D parameter the K factor is shghtly 
above two for the whole range, while for acoplanarity it is 
even larger and increases for larger values of A. This sug- 
gests that A cannot be reliable calculated in perturbation 
theory. 

Finally, in Fig. 3 we compare the next-to-lcading or- 
der QCD prediction for the D parameter to L3 data ob- 
tained at the Z° peak ||l^ and corrected to hadron level. 
The inclusion of the higher order correction decreases the 
discrepancy between the next-to-leading order QCD pre- 
diction and the data. However, there still remains sig- 
nificant discrepancy. This difference may come in part 
from hadronization effects, and also from the uncalcu- 
lated even higher order contributions. 
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FIG. 2. K factor of the D parameter and acoplanarity. 
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FIG. 3. Comparison of the next-to-leading order QCD pre- 
diction for the D parameter differential distribution, ^ g§ to 
L3 data obtained at the peak and corrected to hadron 
level. The upper edge of the theoretical band is obtained 
with renormalization scale a;^ — 0.1, while the lower edge at 
X,. = 1. 



In this Letter we presented for the first time a next-to- 
leading order calculation of the differential cross section 
of two classic four jet shape variables, the D parame- 
ter and acoplanarity. We gave explicit results for the 
radiative corrections to the leading order cross sections. 
The corrections are large indicating that the uncalculated 
even higher order terms are important. This feature is es- 
pecially dramatic in the case of acomplanarity suggesting 
that this observable cannot be reliable calculated in per- 
turbation theary. We also compared the four-jet rates 
obtained by our program to the results of Dixon and 
Signer |l^ and found agreement. 

These results were produced by a partonic Monte Carlo 
program that can be used for the calculation of QCD ra- 
diative corrections to the differential cross section of any 
kind of four-jet observable in electron-positron annihila- 
tion. 

This research was supported in part by the EEC 
Programme " Human Capital and Mobility" , Network 
" Physics at High Energy Colliders" , contract PECO ER- 
BCIPDCT 94 0613 as well as by the Hungarian Scientific 
Research Fund grant OTKA T-016613 and the Research 
Group in Physics of the HAS, Debrecen. 

Note added: After the completion of this work, the 
helicity amplitudes of the e+e^ , 7* — > qqgg pro- 

cess have been published fl^ , and the agreement with 
the results of Ref. in the 7* channel has been es- 
tablished. Also, the authors of Ref. |jll[ pointed out 
a slight error in our binning procedure in the case of 
the Geneva algortihm. Correcting this error we find 
i?4(y^„t = 0.05) = (8.37 ± 0.12) ■ 10-2 that agrees with 
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